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We present a scheme for robust and efficient projection measurement of a qubit consisting of the 
two magnetic sublevels in the electronic ground state of alkali earth metal ions. The scheme is based 
on two stimulated Raman adiabatic passages (STIRAP) involving four partially coherent laser fields. 
We show how the efficiency depends on experimentally relevant parameters: Rabi frequencies, pulse 
widths, laser linewidths, one- and two-photon detunings, residual laser power, laser polarization and 
ion motion. 
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I. INTRODUCTION 

Quantum computation is a promising technique for ef- 
ficient solving of high complexity problems, which are 
inaccessible by classical algorithms [1| . Currently, efforts 
are made within many fields of physics in order to explore 
the possibility of realizing quantum computing. Notable 
examples are superconducting circuits 0, [H .semicon- 
ductors [H, linear optics with single photons Op, [Hi, c old 
neutral atoms in cavities 01 and lattices [cl l9l. llOl] as 
well as cold, trapped ions [11| . So far most prog ress has 
been made in ion tra p sy stems [l~2l . filj , [LH Il5l . [iH li~7| 
where quantum gates |l2l . Il3j , many qubit entanglement 
[HI, [l|| and quantum error correction has been demon- 
strated [!|. 

The qubit states most successfully implemented with 
trapped ions are two hyperfine components of the 9 Be + 
ground state [T3] and a combination of the metastable 
3D 5/ / 2 state and the 4S!/ 2 ground state in 40 Ca + [l3[. In 
both cases linear Paul traps are used to confine the ions. 
In this work we consider the electron spin in the alkali 
earth metal ion ground state as a qubit. In Fig. Q] the 
electronic ground state is denoted by |1) and the qubit 
basis by ||) and \[). Single qubit operations as well as 
gate operations can be performed by driving stimulated 
Raman transitions between the two qubit states 2~l| 



or 



alternatively using controlled mechanical light forces 
Since solely ground states are involved, the qubit deco- 
herence will be limited only by ambient noise fields and 
ion heating effects, rather than excited state lifetimes. 
Limitations due to ambient magnetic field noise can be 
avoided u sing logical qubits of a decoherence-free sub- 
space [23l . l24j . In order to achieve fault tolerant quan- 
tum computation the detection error rate must be kept 
low and different error correction codes estimates allowed 
error rates between 10~ 5 and 10~ 2 [25l . |26| . However, 
high error rates require a large overhead of qubits to en- 
code the error correction. In this paper we present a 
scheme for potentially efficient qubit projective measure- 



ments via shelving of population of one qubit state in 
the long lived metastable D 5 / 2 state found in the iso- 
topes with no nuclear spin of the alkali earth metal ions 
Ca + , Sr + and Ba + as well as the transition metal ion 
Hg+. 

The qubit shelving is performed via a double Stim- 
ulated Raman Adiabatic Passage (STIRAP) process 
US HI] , as illustrated in Fig. [TJ Initial and final states 
of the shelving process are denoted ||) and |5) respec- 
tively. After shelving, the atomic population remaining 
in the 1 1) state can be observed directly by resonantly 
driving the |1) — > |2) and 13) — > |2) transitions and mon- 
itoring the fluorescence [231 . We use the first STIRAP 
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FIG. 1: (Color) Relevant energy levels of isotopes of Ca + , Sr + , 
Ba + and Hg + with no nuclear spin: |1) = IS1/2), |2) = IP1/2), 
|3} = |D 3/2 >, |4) = |P 3/2 ) and |5> = |D 5/2 ). We consider 
the qubit encoded in the Zeeman sublevels, | \) and | J.) of 
the electronic ground state, |1). The states |3) and |5) are 
metastable with lifetimes varying from a fraction of a second 
to many seconds depending on the ion species. The decay of 
these states is neglected here. Since the scheme consists of 
two stages of STIRAP, we break down the five level system 
into two three level A systems, each coupled by two optical 
fields. 



2 



process to transfer the population from ||) to |3), where 
we achieve the spin-state selectivity by using circularly 
right handed polarized light to couple the states ||) and 
|2). The second process transfers the population from 
|3) to the non-fluorescent state, |5). A weak externally 
applied magnetic field defines the quantization axis and 
creates a Zeeman splitting of ||) and |f) (See Fig. [I]). 

STIRAP has been shown to be a robust way of adia- 
batically transferring population from one quantum state 
to another in a three level lambda-system using two laser 
pulses in a counterintuitive order [28|, [3(| HH, HH, HH, H3, 
I35L I36L l37l | . Compared to population transfer via 7r-pulses 
or rapid adiabatic passage, STIRAP has the advantage 
that no strict control of laser amplitude and phase is 
required to maintain a high efficiency [38| . However, 
STIRAP does require the pairs of laser pulses involved 
to be phase coherent relative to each other. In prac- 
tice, this means that either the lasers must be mutually 
phase locked or the laser pulses sufficiently short, such 
that laser decoherence is negligible. The latter method is 
experimentally attractive because high obtainable Rabi 
frequencies permit short pulses to be used and hence 
laser phase locking is circumventable. This approach 
has previously been used in the experimental realizations 
of STIRAP [2!, [37|. Also note that the two STIRAP 
stages can be decoupled within the lifetime of intermedi- 
ate metastable state, 1 3). Of course decoherence is always 
at play. Hence, in the numerical calculations below we 
assume finite laser linewidths and only partially coherent 
lasers. 

The goal of this work is to theoretically investigate 
the shelving process and identify the parameters of im- 
portance to its efficiency. In Sec. |H] we first make some 
general remarks on the STIRAP-process. This analyti- 
cal treatment is valid for any three level A system and it 
serves to develop the criteria for maintaining adiabaticity 
using Gaussian shaped pulses. We then turn to numer- 
ical simulations of the density matrix evolution for the 
full five level system. Here, we consider the roles of one- 
and two-photon dctunings of the Raman resonances, laser 
linewidths, laser pulse widths, Rabi frequencies, stray 
laser light, laser polarizations and ion motion. The re- 
sults are presented in See's ITVl and fVl Here we use as 
an example the 40 Ca + ion, but, apart from different de- 
cay rates and wavelengths, our treatment is also valid for 
Sr+, Ba + and Hg + ions with no nuclear spin. In Sec. [VI] 
we present simulations taking all effects into account and 
finally we conclude in Sec. IVIII 



II. STIRAP WITH GAUSSIAN PULSES 

Let us first consider one STIRAP process involving the 
three atomic basis states (|1), |2), |3)). In this section we 
will neglect spontaneous decay of the |2) state, since this 
is of no importance to the adiabaticity underlying STI- 
RAP. Two mutually coherent and monochromatic laser 
fields A (pump field) and B (Stokes field) couples |I) and 



1 3) to 1 2) (See Fig. [2a)). The interaction Hamiltonian for 
this system in the rotating wave approximation is given 
by 

h ( n A (t) \ 
Hj(t) = - n A (t) 2A A n B (t) (i) 

1 \ n B (t) 2(A A - A B ) J 

where and Ab are the detunings of the two applied 
fields having the real, time dependent Rabi frequencies, 
£lA(t) and ^s(t) respectively. These parameters are de- 
fined as in [2Sj . On two-photon resonance, Aa — Ab = 0, 
diagonalizing |T]) yields the eigenvalues 

w ± = i(A A ±y / A| + n\ + n|) , c/ = o, (2) 

corresponding to eigcnstates 

|+) = sin(a)sin( / 8)|l)+cosC9)|2)-cos(a)sin(/3)|3), 
|-) = sin(a)cos(/3)|l) -sin(/3)|2) + cos(a) cos(/3)|3), 
\d) = cos(a)|l) — sin(a)|3), (3) 

with tan(a) = and tan(/3) =y / -^ T . The |d)-state 

has zero interaction energy and hence it is decoupled 
from the light. As a result adiabatic following of this 
state does not populate of the potentially short-lived |2)- 
state. Now, with all population initially in |1) and only 
the Stokes field applied, this initial state is \d). Adiabat- 
ically decreasing Qb while increasing Qa corresponds to 
a variation of tan(a) from to oo. The pulse sequence 
thus changes the dressed state \d), from |1) to |3). 

Since STIRAP requires adiabatic following of atoms 
while remaining in the |d)-statc, we will now consider 
the criterion for this to happen. We assume essentially 
Fourier limited optical pulses with a Gaussian time de- 
pendence. These give rise to the time varying Rabi fre- 
quencies 




FIG. 2: (a) Three level lambda-system with a pump laser 
(A) and a Stokes laser (B) applied. The two laser fields have 
Rabi frequencies £1a and f2s and they are detuned A a and 
As from resonance, (b) STIRAP pulse sequence. At is the 
pulse separation and r is the full width at 1/e height for both 
pulses. 
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where the peak Rabi frequencies are CIa, o and Qb,q, At 
is the pulse separation and r is the full width at 1/e 
height, assumed to be the same for both pulses (See pulse 
sequence in Fig. H^b)). We parameterize the problem 
in terms of the scaled time, 9 ~ VSt/r, scaled pulse 
separation, 77 = v^At/r and Rabi frequency asymmetry 
r = Ha.o/^b.o- Hence 



= Oi,oexp [-1/2(0-7?) 
n B (6) = fts.oexp [-1/2(0 + rj) 



and 



\d) 



,1/2 779 



|3) 



yj^^>2">i0 _|_ ^— lg — 2r/0 



(5) 



(6) 



Adiabaticity requires the rate of change of the wave func- 
tion to be small compared to the energy separation be- 
tween the dressed state eigenvalues, 



dt 



(7) 



The rate of change of the wave function © is parame- 
terized by the Bloch sphere polar angle in the {|1), |3)} 
basis, which from and ^ is found to be a = 
arctan[rcxp(2770)]. The wave function rate of change is 
then 



da 

~d~9 



2 V 



+ r 



,-l p -2?70 ' 



(8) 



This should be compared to the dimcnsionlcss eigenvalue 
of the energetically closest dressed state, |— ). On two- 
photon resonance and in the limit of large single photon 
detuning, flAfl^B.o "C |Aa|, the eigenvalue is given by 



A_ = 



V8 



= — 2A I 



= 2 " e +r- 1 e- 2 " e )exp [- (9 2 +if)] , 

(9) 

We can now formulate an adiabaticity criterion, which 
should be fulfilled at all times during the pulse sequence, 
where the Rabi frequencies are appreciable. The criterion 
states that 



A(e,r),r) 



I 77 I exp [9 2 



(re 



2tj0 



„-l„-2?70 



< A, 



(10) 



where we identify the parameter relevant for maintaining 
adiabaticity 



A 



16V2|A A | 



(11) 



On inspection A is simply proportional to the product of 
the maximum achievable Raman Rabi frequency and the 
duration of the pulse sequence. Hence, Eq. (fit))) simply 
states that many Raman Rabi cycles should take place 
during the time span of the STIRAP process, analogous 
to the result derived in 12811. 
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FIG. 3: (Color) Maximum value of the adiabaticity function, 
Amax{ri,r) = A(8 max , 77, r) as function of STIRAP pulse sep- 
aration, 7i, for r = 1, r = 2 and r = 4. For r > 1 a solution 
does not exist for smaller rj values, hence the truncation of 
the corresponding curves. The maximum is found by solving 
Eq. (Til}. 



The timing of the STIRAP sequence is contained 
in the function A(9,r],r). For a range of parameters 
this has a maximum as function of time, A max (r],r) = 
A(9 m ax 1 ri,r), where 9 max solves the equation 



2 tanh(2776> ma2; + lnr) = 9 max /r). 



(12) 



For an efficient STIRAP process to occur, given a certain 
value of r, we obviously want to choose the pulse sepa- 
ration, 77, so that A max is small and obeys the inequality 

(ffoD- 

In Fig. [3] we plot A max as function of 77 for values, 
where Eq. (| 1 2[) has a solution. The three curves cor- 
respond to r ~ {1,2,4} and it should be noted that 
AmaxiVi r ) = A max (r),r~ l ). From the graphs we infer 
that values of A in excess of 3-5 should be obtained in 
order to ensure efficient population transfer. Adiabaticity 
is favored when the pulse separation is decreasing toward 
7/ = 0.7 and when the Rabi frequencies are balanced, 
correspondin g to r approaching unity. In an experimen- 
tal situation [37| residual light illuminating the ions be- 
fore and after the STIRAP pulse sequence can become 
a severe limitation to the population transfer efficiency. 
This will be discussed in more detail in Sec. IIVE1 where 
one conclusion is that a short duration of the experi- 
ment is needed in order to avoid population rcpumping 
due to residual light driving a parasitic Raman resonance. 
Hence, we introduce the time duration of the experiment, 
O, which when inserted in Eq. ^ provide the population 
transfer efficiency: 



1 



" 2 exp(-47j6)]" 



(13) 



Obviously, (fl~3"|) shows that the transfer efficiency grows 
with increasing 77, due to the tails of the Gaussian pulses. 
However, for larger 7/ fulfilling the adiabaticity criteria 
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(fTO|) becomes difficult and we hence expect to find an op- 
timum value of the pulse separation, r] opt , which depends 
only on r and A. With respect to r, ^(0) grows for 
high values, while adiabaticity requires values close to 
unity and therefore one might expect an optimal value 
of r above 1 with the exact value depending on A. To 
optimize the transfer efficiency with respect to r\ and r, 
we solve the optical Bloch equations for the three level A 
system. Details of our derivation of the Bloch equations 
can be found in the next section. 

In order to make our Bloch equation solutions general, 
we again ignore the decay from the short lived state |2) 
and we assume perfectly coherent laser fields. The Bloch 
equations are integrated numerically for one-photon de- 
tunings of 4 GHz. These detunings are chosen sufficiently 
large, so that for Rabi frequencies in the 100 MHz range, 
the results are independent of the specific values of the 
detunings. It is checked numerically that increasing the 
detunings even more has no impact on the solutions. In 
this limit our results can be considered general and sys- 
tem independent. 

For fixed values of r and A the population transfer is 
computed as a function of 77, and the optimum value, 
Vopt, corresponding to the maximum population trans- 
fer, P m ax, is found. Examples of such curves can be 
found in Sec. IIV Al where we analyze the sensitivity of 
the population transfer with respect to fluctuations of 77. 
Fig. 01a) shows rj opt as function of r for A = {1,3,10}. 
From the curves it is found that for values of r far from 
unity, a larger pulse separation is preferable. As men- 
tioned above, this is attributed to the wings of the larger 
amplitude Gaussian pulse, which must decay (r < 1) - 
or grow (r > 1) sufficiently before the pulse sequence is 
terminated. This is also the reason why rj opt grows from 
about 0.7 for A = 1 to about 1.1 for A = 10. Fig. g^b) 
shows the population transfer inefficiency, 1 — P ma x, as- 
sociated with the values of rj opt computed in Fig. Ufa) . 
From these curves it is obvious that balanced Rabi fre- 
quencies are crucial for an efficient population transfer 
when the adiabaticity criterion (|10[) is only marginally 
fulfilled, as it is the case for A = 1. As A grows, we see 
that higher asymmetries in the Rabi frequencies can be 
tolerated. This is in accordance with the discussion of 
Fig. [3] It should be noted that the inclusion of finite 
laser lincwidths tends to make rj opt smaller as a result 
of decoherence mechanisms explained in Sec. IIV CI For 
the remaining simulations in this paper we have sought 
to use optimum values of r\ according to the reasoning 
above. 

Although Fig. BJb) shows that the fidelity of qubit 
detection in general can be high for A = 10, it should 
be noted that the difference between a detection error of 
10~ 3 and 10 -5 becomes very important when error cor- 
rection protocols are considered [l[. Hence, optimizing 
the fidelity with respect to r, according to Fig. Q|b), 
could be crucial if the qubit detection scheme should be 
utilized for scalable quantum computation. 




r 

FIG. 4: (Color) (a) Optimum pulse separation, rj op t as func- 
tion of Rabi frequency ratio, r = Qa,o/^b,q- As indicated, 
the three curves correspond to different values of A, given by 
A = 1 (♦), A = 3 (a) and A = 10 (•). One-photon detunings 
are Aa/27t = As/27r = 4 GHz and all results are obtained 
with r — 2 fis. (b) The transfer inefficiency 1 — P ma x cor- 
responding to the values of r\ ovt calculated in (a). The oscil- 
lations of the curves are attributed to Rabi dynamics and it 
should be noted that the curves are symmetric with respect 
to letting r — > r~ . 



III. OPTICAL BLOCH EQUATIONS 

The previous section developed analytical criteria for 
maintaining adiabaticity when Gaussian pulses are used 
and when the one-photon detunings are much larger than 
Rabi frequencies. We now turn to numerical simulations 
of the detection scheme taking into account the full dy- 
namics and using realistic experimental parameters. We 
solve the optical Bloch equations describing the dynamics 
of the five levels depicted in Fig. [TJ but ignore magnetic 
sublevels. We stay within the dipole approximation and 
define the Rabi frequencies 

n A = Li 12 E A /h (14) 
fis = mE B /n (15) 
n c = mEo/h (16) 

&d = fJ-5iE D /h (17) 
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where faj are the dipole matrix elements and Ei laser 
field amplitudes. The laser fields with frequencies Wij 
are not necessarily on resonance and we introduce de- 
tunings A.,j = ' h 3 — uiij with respect to the atomic 
energy levels {Si}. With the notation of Fig. Q] we define 
A A = A 2 i,A s = A 23 ,Ac = A 43 and Ad = A45. Fur- 
thermore, we apply the rotating wave approximation to 
arrive at the interaction Hamiltonian for the system 

5 

H = ^ hLU i p li + (hQ A Pl2 + hD.Bp32 + fi^CP34 + flQ D p 5 4 + h.C.), 
i=l 

(18) 

where u>i = £i/fr and the density operator elements are 
Pi] = b)(*l- Spontaneous emission is introduced as de- 
cay terms in the density matrix elements. We neglect 
spontaneous decay from the metastable |3) and |5) states 
because their lifetimes arc significantly longer than the 
simulation time (See lifetimes of the different ion species 
in appendix). Phase fluctuations of the lasers are intro- 
duced as decay of the coherences and ion micromotion as 
a harmonic modulation of the detunings. 

In the simulations we use the specific wavelengths 
(X A = 397 nm, \ B = 866 nm, \ c = 850 nm, Ad =854 
nm) and spontaneous decay rates {Y~2i/^ = 21 MHz, 
r 23 /27r = 1.7 MHz, T ai /2k = 22 MHz, r 43 /27r = 0.18 
MHz, r 45 /27r = 1.6 MHz) for the 40 Ca+ ion, but this 
apart the calculations are also valid for isotopes of Sr + , 
Ba + and Hg + with no nuclear spin (See relevant pa- 
rameters for the most abundant isotopes in appendix). 
We use Gaussian pulses as defined in Eq. (j4|) with pulse 
widths Ti = 2/2S. All other parameters are varied depend- 
ing on which investigations are made. When there is no 
argument for a different choice we use Rabi frequencies 
fl h0 /2TT = 100 MHz and detunings A ? ;/2tt = 600 MHz. 
The delay between pulses is chosen optimal, which in the 
simulations presented is either At = 1.2/xs or At = l.3[is. 
The theory of Sec. [H] applies when Ifi^ol |Aj|, which 
is the case in the majority of the simulations. In these 
cases we specify the values of ry and A to make compar- 
ison with Sec. |TT] straightforward. For the parameters 
mentioned above 77 = 0.85 or 0.92 and A = 9.3. All pa- 
rameters for the simulations are mentioned in the figure 
captions and hence omitted in the text. 

To recapitulate, the detection scheme consists of two 
STIRAP processes. The first transfers population from 
|1) to |3) applying the A and B-fields as pump and Stokes 
field, respectively. The second process between |3) and 
1 5) uses C as the pump field and D as the Stokes field. 

We first investigate the effect of pulse delay (Sec. lIV All , 
laser detunings (Sec. HVBj) . laser linewidth (Sec. IIV Cf . 
pulse width (Sec. lTVD|) and residual light (Sec. lTVE|) . All 
these simulations are performed on the latter transition 
because this transition is subject to spontaneous decay 
out of the sub-system and thus has a higher sensitivity to 
non-adiabaticity. An experimental investigation of this 
transition can be found in [37| . 

The first STIRAP stage is more sensitive to ion mo- 
tion because the two laser fields used have very different 
wavelengths and the Dopplcr shift induced by the ion mo- 



tion therefore gives a large two-photon detuning. Since 
this stage is responsible for the internal state selection, 
it is also important to investigate the effect of laser po- 
larization errors and the resulting depletion of the wrong 
qubit state. The simulations of ion micromotion (Sec. 
HVFp and polarization errors (Sec. [V]) are therefore per- 
formed on the first STIRAP transition. The simulations 
are made with all the population initially in |1) or |3) 
and we find the transfer efficiency as the final population 
in 1 3) or |5) for the first or the second STIRAP stage 
respectively. We assume that the calculated transfer effi- 
ciency is equivalent to the qubit detection efficiency, and 
we do thereby not focus on signal to noise issues of the 
fluorescence detection of the non-shelved state ||). 



IV. SIMULATIONS IGNORING ATOMIC SPIN 

In this chapter we assume perfect laser polarization 
and hence ignore irrelevant magnetic sublevels of the 
ions. As mentioned above, we start out by considering 
the latter |3) |4) -> |5) STIRAP process. 



A. Effect of pulse delay 

The first parameter we investigate is the delay between 
the two pulses and in Fig. [5] we show the transfer effi- 
ciency as a function of -q for various peak Rabi frequen- 
cies Vtcfi — Qd.o- Both laser fields are on resonance. 
The simulations show that for Rabi frequencies above 




FIG. 5: (Color) Transfer efficiency as a function of delay be- 
tween pulses for different choices of peak Rabi frequencies: 

( ) Slc.o = tt D ,o = 2tt x 10 MHz, (■■■) flc.o = Od, = 

2tt x 20 MHz, ( ) n c ,o = Od,o = 2tt x 100 MHz, (— ) 

Qc,o = fio.o = 2-7T x 300 MHz. Positive delay correspond to 
the counter intuitive pulse sequence, where the laser pulse of 
field D arrives before the laser pulse of field C. Parameters: 
tc — td = 2 (is and Ac = Ad =0. 

ftc.o = ^d.o = 27T x 100 MHz the transfer efficiency is 
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close to unity and insensitive to fluctuations in the delay 
(P 5 > 0.999 for fluctuations of ij less than 0.2 around 
the optimal value of 0.85). For smaller Rabi frequencies 
the transfer efficiency decreases and the sensitivity to de- 
lay fluctuations increases. The plateau of P5 ~ 0.07 for 
positive r] values is the result of optical pumping by the 
pump field. Qc- For r\ > this pulse is applied last, and 
the Stokes field cannot repump population. 

We wish to investigate how changes in the delay be- 
tween pulses influence the transfer efficiency when it is 
not diminished by other effects. The results of Fig. 
are therefore calculated for zero one-photon detuning. 
In this limit the theory of Sec. [TT| is not valid, but in 
the rest of the work one-photon detunings considerably 
larger than Rabi frequencies will be used and hence the 
analytical theory can be used to interpret the numerical 
simulations. 

The results presented in Fig. [5] still show the same 
trends as the analytical theory. High Rabi frequencies 
are needed to maintain adiabaticity and, as indicated by 
the crosses, the delay giving maximum transfer efficiency 
is found to grow with increasing Rabi frequencies, con- 
sistent with Fig. [U 



B. Laser detunings 

To reduce the probability of incoherent diabatic exci- 
tations we introduce large one-photon detunings of the 
lasers. Thus, we now investigate the sensitivity of STF 
RAP to one- and two-photon detunings. The simula- 
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FIG. 6: (Color) Transfer efficiency as a function of two-photon 
detuning (Ad-Ac) of the C and D lasers for different choices 
of one-photon detuning, Ac- The simulations are made with 
fic.o = &d,o = 27r x 100 MHz, tq = td = 2 fis and At = 1.2 

Ms (77 = 0.85). ( ) Ac = MHz, ( ) Ac = 2tt x 300 

MHz (A = 18.5), (■■■■) Ac = 2tt x 600 MHz (A = 9.3), (-•-•) 
Ac = 2tt x 1200 MHz (A = 4.6). 

tions, presented in Fig. [SJ show that increasing the one- 
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FIG. 7: (Color) (a) Transfer efficiency as a function of two- 
photon detuning. Evolution of eigenvalues as a function of 
time is shown in (b):A D -Ac = MHz, (c) Ad-Ac = 27rx2 
MHz and (d) Ad — Ac = -2tt x 2 MHz . In (b), (c) and (d) 

we use the signatures: ( ) uj d , (- - -) uj~ . Parameters 

used for all graphs: Q c ,o = 2-k x 100 MHz, Q n ,o = 2tt x 200 
MHz, Ac = 2tt x 600 MHz, rc = td = 2 /is and At = 1.2 /Lis. 
(A= 18.5, 77 = 0.85). 



photon detuning does not limit the transfer efficiency 
as long as we are close to the two-photon resonance 
An — Ac = 0. This criterion gets stricter as we increase 
the one-photon detuning. If we require P5 > 0.99 the de- 
mand on two-photon detuning is Ad — Ac± 2tt x 0.5 
MHz when Ac = 2n x 1200 MHz as estimated from 
the dash-dotted curve of Fig. El For Ac = 2tt x 600 
MHz, maintaining a two-photon resonance within 2tt x 1 
MHz is sufficient as seen from the dotted curve of Fig. [5] 
An efficient transfer can thus be maintained in spite of 
a drift from two-photon resonance making STIRAP ro- 
bust with respect to small frequency drifts of the involved 
lasers. A closer look at Fig. [6] reveals that the spec- 
tra are not exactly symmetric with respect to the sign 
of the two-photon detuning. With unbalanced Rabi fre- 
quencies, ttD,o = 2^c.o (See. Fig. 0Ja)), the asymmetry 
becomes more evident. For small negative two-photon 
detunings the transfer efficiency is much higher than for 
the corresponding positive two-photon detuning. This ef- 
fect is due to diabatic transfer between \d) and the ener- 
getically closest bright state |— ). As discussed in Sec. |TT] 
this is likely when \d) and |— ) are nearly degenerate. The 
eigenvalues found in Eq. @ are valid only on two-photon 
resonance, and the solutions are more complicated in the 
situation, where this condition is not met. A general ex- 
pression for the eigenvalues has been derived in [39| and 
is plotted in Fig. [7{b-d). Here, we show the eigenvalues 
of \d) (solid curves) and |— ) (dashed curves) as a function 
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of time for three different choices of two-photon detun- 
ing. In the case where A D - A c = MHz (See Fig. Eft))) 
the eigenvalues of course coincide when the Rabi frequen- 
cies are zero before and after the pulses, but in this case 
no diabatic transfer will occur as no coupling is present. 
For positive and negative detunings (Ad — Ac = 2ir x 2 
MHz, Fig. E[c) and Ad — Ac = —2ir x 2 MHz, Fig. Eld)) 
we see avoided crossings leading to a probability for di- 
abatic transfer to the |— )-state. Such a transfer leads 
to population of the |4)-state, which decays rapidly, and 
the population mainly ends up in |1). When the energy 
splitting between the \d) and the |— ) state is small, main- 
taining adiabaticity requires high coupling strengths. For 
a positive two-photon detuning the avoided crossing oc- 
curs late in the sequence, as shown in Fig. Elc). Here 
the coupling strength is relatively small due to the asym- 
metry Qd,o — 2f2c.0 7 hence the probability of diabatic 
transfer is large. With a negative two-photon detuning 
the avoided crossing occurs early, as shown in Fig. Elc). 
The coupling is stronger at this time and the probabil- 
ity of diabatic transfer is thus smaller. This difference 
gives rise to the asymmetry of the two-photon spectrum 
in Fig.Ela). 
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FIG. 8: (Color) Transfer efficiency as a function of the 
linewidth of the C and D laser fields for different choices 
of Rabi frequencies. The simulations use tc = td = 2 fis, 

At = 1.2 ^s (?7 = 0.85) and A c = A D = 2tt x 600 MHz. ( ) 

P 5 = l/3 + 2/3exp(-3735T 2 /i6At), ( )n co = Q D = 

2tt x 300 MHz (A = 83.3), (■■■) n c ,o = ^d,o = 2tt x 100 MHz 
(A = 9.3). 



C. Laser linewidths 

Laser phase fluctuations lead to dephasing between the 
three states involved and may therefore affect STIRAP. 
This dephasing effect has been studied thoroughly in [|(| , 
where the STIRAP transfer efficiency has been shown to 
depend only on the dephasing rate between |3) and |5), 

735- 

p 5= l + 2 e -3 735 r 2 /16A t; (lg) 

for Gaussian pulses within the adiabatic limit, not taking 
decay from |4) into account. In Fig. [5] we present sim- 
ulations showing the exact influence of the laser phase 
fluctuations parameterized by the laser linewidth. Also 
shown as the full curve is Eq. (fT9| . The curves show 
that the transfer efficiency is strongly limited by the 
laser phase fluctuations and because of the non-zero de- 
cay from 1 4) the situation is even worse than predicted 
by Eq. (fl9)) . For high Rabi frequencies the coupling is 
stronger and therefore the limitations due to dephasing 
are less pronounced. To maintain a transfer efficiency 
above 0.99 we find laser linewidths required to be below 
2tt x 1.5 kHz for Q c ,o = ^d,o = 2n x 100 MHz and below 
2?r x 2 kHz for Q c ,o = ^d,o = 2ir x 300 MHz. It should 
be noted that by laser linewidth we mean the frequency 
fluctuations averaged over the time duration of the ex- 
periment. This is not necessarily the steady state laser 
linewidth. For pulse sequences of roughly 10 /is duration, 
as are employed in our simulations, laser linewidths of a 
few kHz are not unrealistic for most laser systems. 



D. Role of pulse widths 

As discussed above, the laser phase fluctuations limit 
the transfer efficiency, making it preferable to use short 
pulses, since this typically reduces the effective laser 
linewidth. Short pulses, however, require higher Rabi fre- 
quencies to maintain the adiabaticity. In Fig. [HI wc show 
the transfer efficiency as a function of pulse width in the 
case of no laser linewidth as well as the case of a 2ir x 2 
kHz linewidth. As expected, for the case of no linewidth 
we find an increasing efficiency as we increase the pulse 
widths because the criterion for adiabaticity is better ful- 
filled for larger pulse widths as found in (TIT]) . We indicate 
the point where we exceed 0.995 on the graph. When we 
introduce a 27T x 2 kHz laser linewidth an optimum pulse 
width is found as indicated on the graph. 

This optimum depends on the Rabi frequencies of the 
lasers and the dependency is pictured in Fig. [TO] In the 
case of no linewidth, instead of the optimum, wc plot 
the pulse widths required to exceed a 0.995 transfer ef- 
ficiency. This result is compared to the optimum pulse 
width, when lasers have 2ir x 2 kHz linewidth. Both 
curves show that short pulses are preferred as Rabi fre- 
quencies grow. However, the mutual dephasing of the 
lasers involved, results in smaller optimum pulse widths 
when finite laser linewidths are introduced. For higher 
Rabi frequencies, the two curves cross, so that shorter 
pulse widths apparently are preferred for the perfectly 
coherent lasers. This is a result of the artificial limit 0.995 
chosen to define the black curve. We have assumed the 
laser linewidths to be constant over the range of pulse du- 
rations involved. More realistically one would expect the 
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FIG. 9: (Color) Transfer efficiency as a function of pulse width 

of the C and D laser fields for 2n x 2 kHz ( ) as well as 

no laser linewidth (- - -) The delay between the pulses is 
given by r\ = 0.85, laser detunings are Ac = Ad = 2n x 600 
MHz and peak Rabi frequencies f2c,o = fiD.o — 2n x 100 MHz 
(A = 9.3). 
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FIG. 10: (Color) Optimal pulse width as a function of the 
Rabi frequencies of the C and D laser fields for 2n x 2 kHz (•) 
as well as no laser linewidth (■). The delay between the pulses 
is given by rj — 0.85, laser detunings are Ac = Ad = 2ir x 600 
MHz. 



laser linewidths to grow proportional to (tc.d) q , where 
the power, a, is in the range 0.5 — > 1 for a phase diffusion 
process. 



E. Residual light 

In the first generation of STIRAP experiments beams 
of atoms passed through two stationary laser beams 
[28l . [30j . Here, the detection region was spatially well 



separated from the STIRAP region and in this case no 
residual light was present. When we consider the situ- 
ation with stationary ions and pulsed laser beams, the 
situation is less favorable. In the case of a constant back- 
ground level of laser light, the ions do not start out in 
an exact dark state and hence population transfer may 
be limited by the resulting finite population of the short 
lived states \2) and |4). Moreover, since the residual light 
is able to excite the Raman resonance between the states 
1 3) and |5), repumping of the transferred population can 
take place after ended STIRAP pulse sequence. With a 
small background level of light this Raman resonance will 
be narrow and hence the limitation of population transfer 
is most severe near two-photon resonance as experimen- 
tally found in [3lJ • Assuming Rabi frequencies of 2n x 100 
MHz and extinction ratio w 10~ 4 of the optical pulse 
generators, we arrive at residual Rabi frequencies in the 
MHz range, which can easily repump a significant portion 
of the shelved population on two-photon resonance. In 
this section we address this problem and establish upper 
limits on the peak Rabi frequencies involved. 

We denote the Rabi frequency of the residual light 
Q fi,i, i = {A, B, C, D}. We first look at the transfer ef- 
ficiency as a function of the two-photon detuning. When 
no residual light is present the transfer efficiency is op- 
timal on two-photon resonance as shown previously in 
Fig. [51 but this is no longer the case when we introduce 
a fraction of residual light, il o s,i/^i,0- I n Fig. HIT a) we 
present the transfer efficiency as a function of two-photon 
detuning for different fractions of residual light. As we 
increase this fraction we see an increasing loss of pop- 
ulation on two-photon resonance. With f2 o ff,i/^i,0 l ess 
than 0.02 (dashed curves) this loss is only a few percent 
relative to the maximum transfer efficiency. However, al- 
ready with f2ofr,i/^i,o = 0.05 (solid curve) the loss has 
grown to nearly 50 % near two-photon resonance, ren- 
dering the detection scheme useless in this case. 

As mentioned earlier, residual light is typically present 
as a fraction of the peak Rabi frequencies, and hence very 
high Rabi frequencies are not necessarily preferable. In 
Fig. [TTTb) we show the transfer efficiency as a function of 
the peak Rabi frequencies in the STIRAP pulses for dif- 
ferent values of the fraction of residual light. Increasing 
the Rabi frequencies increases the transfer efficiency un- 
til the residual light limits the transfer. The maximum 
transfer efficiency is found to be 0.99 around the Rabi 
frequencies, flc,o = ^d,q — 2tt x 50 MHz with a residual 
light level less than f2 o ff,i/^i,0 = 0.02. For even higher 
Rabi frequencies we observe Rabi-floppings between |3) 
and 1 5). 

Experimentally VLifi/2iv =100 MHz can be achieved 
focusing a beam with a mW power level to a fiin spot size. 
Switching of the lasers with acusto-optical modulators or 
electro-optical modulators can be done on ns timescale 
reducing the power level to a few hundred nW which cor- 
respond to £l o &,i/£li,0 — 0.01. This is represented by the 
green dotted curve in Fig. IllT b) , which obviously reveals 
the limitations set by residual light. The extinction ra- 
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FIG. 11: (Color) (a) Transfer efficiency as a function of two- 
photon detuning with f2c,o = f2zj,o = 2tt x 100 MHz and 
A c = 2tt x 600 MHz. A = 9.3. (b) Transfer efficiency as 
the peak Rabi frequencies are varied, simulated with Ac = 
Ad — 2tt x 600 MHz. The different curves represent different 
fractions of residual light: (•••) fi o ff,i/fii,0 = 0.01, (- - -) 

fJoff,i/f2i,o = 0.02 and ( ) Q o K,i/£li,0 = 0.05. Parameters 

used in simulations: rc = td — 2 fis, At = 1.2 [is, yielding 
r\ — 0.85. The total time duration is 20 /is and laser linewidth 
has been ignored. 

tio 10 -4 correspond to the situation where the lasers are 
extinguished by one optical pulse generator. Applying 
additional AOM's and EOM's in succession will reduce 
the residual light level, and as a result, the limitation 
due to this light can be smaller as well, at the expense of 
increased experimental complexity. 



F. Ion motion 

The ions have so far been considered as having no mo- 
tion with respect to the laser fields. We consider the ions 
to be confined in a linear Paul trap [4lj , where they ex- 
perience an effective confining potential due to RF- and 
DC-ficlds. The ion motion can be derived to consist of a 
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FIG. 12: (Color) Transfer efficiency as a function of the micro- 
motion velocity. The different curves show results for differ- 
ent detunings: (■■■) A A = A B = MHz, ( ) A A = A B = 

2tt x 300 MHz (A = 18.5), ( ) A A = A B = 2tt x 600 MHz 

(A = 9.3) and (— ) A A = A B = 2tt x 1200 MHz (A = 4.6). 
Parameters used for simulations: ta = r B = 2 fis, At = f .3 
/is (77 = 0.92), ft A = Q B = 2ttx100 MHz and = 2ttx16.8 
MHz. The temperature of the ion is defined as T=mv 2 /k B , 
where m is the ion mass and k B Boltzmann's constant. The 
Doppler temperature for Ca + , Td, is indicated by the vertical 
line. 

three-dimensional harmonic motion with frequencies to x , 
u) y and uj z called the secular motion [4lj . In addition 
the RF-field drives an oscillation with the RF-frequency, 
Qrf- This oscillation is called the ion micromotion. The 
secular frequencies are typical much smaller than Qrf 
and hence STIRAP transfer efficiency is mainly limited 
by the micromotion. The effect of micromotion is simu- 
lated by modulating the detunings as 

Ai(t) = A<(0) - v^- cos(n RF t). (20) 

We only consider micromotion in one dimension with 
maximum velocity v. This classical treatment of micro- 
motion is justified by quantum Monte Carlo simulations 
where the ion motion has been quantized showing no vis- 
ible difference between classical and quantized motion 
[42| . Another argument for this classical simplification 
is the strength of the laser fields making fl i ^> FIrf- 
As discussed previously, the STIRAP transfer efficiency 
is more sensitive to two-photon than one-photon detun- 
ings. If we assume co-propagating laser fields and take 
the Ca + ion as an example, the two-photon detuning 
will be strongly modulated in the case of the |1)-|2)-|3) 
transition where the two laser fields have a large wave- 
length difference (AA/A) i23 = 0.54, while the |3)-|4)-|5) 
transition experiences a much smaller modulation of the 
two-photon detuning (AA/A)345 = 0.005. This means 
that the first transition, in the case of 40 Ca + , is strongly 
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limited by micromotion while the second remains virtu- 
ally unaffected. 

In Fig. [12] we show the transfer efficiency as a function 
of the micromotion velocity for different choices of de- 
tuning. The simulations show that when we increase the 
one-photon detuning the micromotion velocity becomes 
an increasingly limiting factor. For small micromotion 
velocities we see that the transfer efficiency is higher for 
small one-photon detunings. This is due to a decreased 
sensitivity to a non-vanishing two-photon detuning as 
discussed in Sec. HVBl For higher velocities the one- 
photon detuning is strongly modulated. This results in 
sidebands, of which the high frequency component comes 
close to resonance, hereby inducing real transitions to the 
short lived |2)-statc. This effect gets more evident when 
the one-photon detuning is close to Vp, which for our 
parameters varies between and 100 MHz. Fig. [T2l shows 
small sensitivity to micromotion for very small and very 
large values of A^ = Ab, while the A a = Ab = 27rx300 
MHz trace shows a rapid drop as a function of v. For 
the chosen parameters a micromotion velocity below 1 
m/s ensures a transfer efficiency above 0.998. From the 
kinetic energy of the ion a temperature can be defined 
through the relation i/c#T= ^mv 2 . Hence, in the case 
of 40 Ca + a velocity below 1 m/s correspond to a tempera- 
ture below 4.8 mK. When the ion is Doppler laser cooled, 
temperatures below the so-called Doppler temperature, 
Td, cannot be obtained. For 40 Ca + To = 0.5 mK as 
indicated by the vertical line on Fig. [T2] This is well 
below the required 4.8 mK and Doppler cooling is thus 
sufficient to achieve a high transfer efficiency. Doppler 
temperatures for other relevant ions can be found in Ta- 
ble D 



V. POLARIZATION REQUIREMENTS 

Until this point Zeeman sublevels have been ignored. 
However, in order to investigate the effect of laser polar- 
ization, deviating from purely circular, we now introduce 
all of the Zeeman sublevels of the first stage STIRAP. 
The optical Bloch equations for the full eight level sys- 
tem are derived analogous to the equations in Sec. IIIII 
We consider the field A to have both a± polarization 
components as well as a possible n component in case of 
magnetic field errors. The second STIRAP field, B, is 
assumed to have perfect <r_ polarization. We denote the 
er_i_ polarized component of VLa by fTj[, the cr_ polarized 
component by fl A and the it polarized component by fl A 
(See Fig. 

Efficient qubit projection requires all population 
moved from \l) to the -D3/2 levels while leaving the |f) 
population unchanged. In order to test these criteria we 
simulate the problem with two different initial states cor- 
responding to all population in |j) and jf) respectively. 
The results are shown in Fig. [TJ] First we focus on the 
situation, where all initial population is in the ||) qubit 
state. Ideally the final populations in both qubit states 




FIG. 13: (Color) Zeeman sublevels of the three states involved 
in the first STIRAP process. The S1/2 and P1/2 states are 
coupled by a laser field with a cr+ polarized component (f2j[), 
a er_ polarized component (fi^ ) and a n polarized component 
(f2^). The D3/2 and P1/2 states are coupled by a laser field 
(£Ib) with a pure <r_ polarization. 



should be vanishing for our scheme to work. The results, 
shown in Fig. [TJ](a) and (b), reveal that detection errors 
can be kept below 0.01 provided the magnetic field can 
be controlled accurately enough to keep Sl^/fij^ below 
0.02. However, when all initial population is in the |t) 
state, the situation is worse. From Fig. [TJ] (c), we find 
that the population in the | f) state is pumped out as 
the relative amplitude of the er_ polarized component is 
increasing. This is of course undesirable. Hence, in order 
to keep qubit detection errors below 0.01 for an arbitrary 
initial state, we must demand that n A /il\ < 0.02 and 
{l A /ft\ < 0.04 with the parameters used in this simula- 
tion. 



VI. FULL SIMULATION 

Finally, we present the results of a simulation of the 
full two stage STIRAP population transfer including all 
of the above described effects except polarization errors. 
The parameters relevant to the simulation can be found 
in the caption of Fig. [151 The figure shows the popu- 
lation transfer efficiency as a function of the two-photon 
detuning, Ab — A a, for the first STIRAP stage. Assum- 
ing that the ions can be kept cooled close to the Doppler 
temperature, we have chosen a temperature of 0.8 mK 
due to the micromotion for this simulation. The result- 
ing peak velocity of 0.4 m/s results in a loss of population 
transfer of less than 0.01 according to the discussion of 
Sec. IIVF1 The two-photon detuning of the second STI- 
RAP stage has been chosen to be zero and the residual 
power level is chosen to yield a Rabi frequency of 0.01 
relative to the peak Rabi frequencies for all four optical 
fields. 

From Fig. [I])] we find, when all effects are taken 
into account, a maximum transfer efficiency of 0.97 for 
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FIG. 14: (Color) Effect of polarization error and error in 
magnetic field direction. Populations of the qubit states 
IJ.), P-, and |t), P+, plotted as a function of the relative 
<j_ polarization component, Q A /Q A . The various curves on 
each plot correspond to different magnetic field direction er- 
rors, parameterized by the relative -n polarization component, 

nyfrt. (— ) n°/n+ = o, ( ) a° A /n+ = o.oi, (- 

— ) fii/fij; = 0.03 and (•••) 0.° A /Q+ = 0.05. The upper 
two graphs, (a) and (b), correspond to all initial population 
in the | J.) qubit state, while the lower two graphs, (c) and 
(d), correspond to all initial population in The param- 
eters for the simulations are 

n A,o = ^b.o = 2tt x 300 MHz, 
A A = A fl = 2tt x 300 MHz, r = 2 /is and At = 1.3 fxs. 
(A = 167, n = 0.92). 
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FIG. 15: (Color) Full simulation of the two stage STIRAP. 
The plot shows the population transfer as a function of two- 
photon detuning of the first STIRAP stage. Parameters are: 
A A = Ac = A D = 2tt x 600 MHz, fl A = Q B = Q c = 
tt D = 27T x 100 MHz, t a = t b = t c = r D = 2/is. All laser 
linewidths are 2n x 2 kHz and the trap is characterized by 
£Irf — 2tt x 16.7 MHz and v = 0.4 m/s, corresponding to 
a temperature of 0.8 mK. The STIRAP's have At = 1.2/is 
corresponding to r\ — 0.85 and A = 9.3 for both stages. The 
two STIRAP stages are separated by 10 fj,s and the overall 
simulation time is 30 ^ts. The Rabi frequencies associated 
with the residual power are 0.01 relative to the peak Rabi 

frequencies for all four fields. The curves correspond to: ( 

) all effects on, (- - -) residual power set to zero but laser 
linewidth on, and (•••) both laser linewidth and residual power 
set to zero. 



Ab-Aa = -2ttx0.09MHz. This corresponds to the full, 
black curve. In order to identify the processes limiting 
this efficiency we show with the green, dashed curve the 
resulting efficiency when the residual optical power is set 
to zero. Comparison with the full curve shows that a loss 
in efficiency of almost 0.01 can be attributed to repump- 
ing by residual laser light. The blue, dotted curve of Fig. 
[Pol shows the transfer efficiency when both residual light 
and laser linewidths are set to zero. Here a maximum 
transfer efficiency of 0.996 is found near two-photon res- 
onance indicating that the finite laser linewidth is indeed 
a deleterious effect to STIRAP. The linewidth of 2 kHz 
assigned to all lasers in the simulation results in a loss of 
transfer efficiency of nearly 0.02. In order to reduce the 
effect of laser phase fluctuations shorter pulse durations 
and consequently higher Rabi frequencies must be used. 

The simulation illustrates the importance of efficient 
switching on and off laser power in the STIRAP process. 
Unless the lasers involved in the two STIRAP stages are 
pairwise phase locked one obviously needs fairly high in- 
tensities of the optical pulses and hence the requirements 
to the residual laser power become more stringent. How- 
ever, as can be seen from Fig. [SJ it is noted that the 
scheme is insensitive to fluctuations in the laser power 
as long as the transitions involved are well saturated and 
adiabaticity is maintained. Experimentally, high Rabi 
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frequencies can be achieved at the wavelength listed in 
the appendix, and the complication of phase locking can 
thus be avoided [37l |. 

VII. SUMMARY AND CONCLUSIONS 

In summary, we analyze the scheme where two mag- 
netic sublevcls of the S1/2 ground state of alkali earth 
metal ions arc considered as qubit states. We have shown 
that qubit projection can be effectively performed by 
electronically shelving the population of one magnetic 
sublevel via a double STIRAP process. Hereby the pop- 
ulation is transferred to the metastable D 5 / 2 state and 
remaining population of the other qubit state can be de- 
tected by driving the strong Doppler cooling transition 
of the ions and spatially monitoring the fluorescence. In 
order to selectively shelve only one qubit state we apply 
appropriately polarized light in the first STIRAP process. 

We have established the formalism describing the sys- 
tem consisting of the 5 energetically lowest levels of the 
alkali earth metal ions and we model the STIRAP pro- 
cesses using Gaussian shaped pulses from partially phase 
coherent lasers. Specific simulations have been performed 
by solving the optical Bloch equations for the 40 Ca + ion. 
From these simulations we conclude that our scheme is in- 
deed feasible with a projection efficiency exceeding 99 %, 
provided laser linewidths can be kept below 1.5 kHz over 
the timescale of the pulse sequence. The short timespan 
required to avoid decoherence of the lasers makes it nec- 
essary to have high Rabi frequencies in the 100's of MHz 
regime in order to maintain adiabaticity. However, since 
we are driving strong, dipole allowed transitions, such 
high Rabi frequencies can be achieved by using moder- 
ately focused lasers with 5-10 mW power. 

In order to avoid unwanted excitations, we assume all 
lasers to be detuned roughly 0.5-1 GHz from the one- 
photon resonance. In the limit of this detuning domi- 
nating over Rabi frequencies, we have identified the adi- 
abaticity criteria for the ions when Gaussian light pulses 
are applied. Given an available laser power and hence 
certain Rabi frequencies, we find an optimum pulse sep- 
aration depending only on the parameters, A and r. 
Within variations of 20 % we find the optimum pulse 
separation to be where r is the l/e full width of 

the Gaussian pulses. 

In contrast to early STIRAP experiments, where 
atomic beams and stationary laser beams were used, our 
system involves stationary atoms and laser pulses. Since 
the detection zone and the STIRAP zone are spatially 
identical for stationary ions, we have to take into ac- 
count the finite extinction ratio of the optical pulse gen- 
erators involved in the experiment. Due to a finite ex- 
tinction ratio, a fraction of the peak Rabi frequencies of 
the STIRAP pulses will be present before and after the 
pulse sequence. As a result, shelved population can be 



repumped, and consequently the shelving efficiency re- 
duced. Clearly mechanical shutters will efficiently block 
laser beams, exposing the ions, but unfortunately such 
shutters need several tens of microseconds to activate. 
On the shorter time scale AOM's and EOM's can be em- 
ployed with extinction ratios around 1 : 10 4 . Hence, we 
have analyzed the role of residual light at the level of 
10~ 2 in Rabi frequency, and this was found to restrict 
the shelving efficiency severely if the Rabi frequency ex- 
ceeded 2ir x 100 MHz. In order for the detection scheme 
to be more efficient, efforts must be made to extinguish 
the light even further. This could be achieved by intro- 
ducing more AOM's or EOM's in succession. 

Motion in the ion trap is included in the analysis, and 
found to be unimportant close to the Doppler cooling 
limit. For calcium ions this is 0.5 mK. 

Finally, we have analyzed the effect of impure polariza- 
tion of light, which could arise from stray birefringence or 
variations in the direction of an external magnetic field. 
Solutions to the optical Bloch equations involving mag- 
netic sublevels of the 3 lowest states of the calcium ion 
show that projection errors can be kept below 0.01 pro- 
vided the relative Rabi frequencies of the unwanted po- 
larization components can be kept below 0.02-0.04. 

A simulation including all effects outlined above, ex- 
cept polarization errors, show the importance of being 
able to efficiently switch on and off laser power as well 
as using highly coherent or perhaps phase locked lasers 
in this qubit detection scheme. Since we assume that the 
lasers involved in the scheme are only partially coher- 
ent, a loss of transfer efficiency of almost 0.02 is found 
for laser linewidths of 2 kHz over the time of the experi- 
ment. Performing the shelving on a shorter timescale will 
typically reduce the relevant linewidths but also require 
higher Rabi frequencies and consequently better suppres- 
sion of residual laser light. 



APPENDIX: RELEVANT ION DATA 

Table[T]shows spectroscopic data for some selected ions 
possessing the level structure and transitions shown in 
Fig. [U Listed are wavelengths, Xi, associated with the 
fields fli, i = {A, B, C, D}, decay rates of the short lived 
states 1 2) and |4) into the stable and metastable states 
as well as decay rates of the metastable states into the 
stable states and the Doppler temperature, Tjj. The data 
is reproduced from Ref . [43[ . 
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TABLE I: Data of relevant ions 



Property 


40 Ca+ 


88 Sr + 


138 Ba+ 


202 Hg+ 


\ A [nm] 


397 


422 


493 


194 


As [nm] 


866 


1092 


650 


10670 


Ac [nm] 


850 


1004 


585 


991 


\ D [nm] 


854 


1033 


614 


398 


r 2 i/27r [MHz] 


21 


20 


14 


69 


r 23 /27r [MHz] 


1.7 


1.5 


5.3 


0.05 


r 4 i/27r [MHz] 


22 


23 


19 


168 


r 43 /27r [MHz] 


0.18 


0.18 


0.76 


0.48 


r 45 /27r [MHz] 


1.6 


1.4 


5.9 


40 


r 3 i/27r [Hz] 


0.15 


0.40 


0.009 


1.62 


r 5 i/27r [Hz] 


0.15 


0.46 


0.003 


7.96 


T D [mK] 


0.50 


0.49 


0.35 


1.7 
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